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Abstract 



In this paper, the influence of aggregate size and volume fraction on shrinkage induced micro-cracking and 
permeability of concrete and mortar was investigated. Nonlinear finite element analyses of model concrete 
and mortar specimens with regular and random aggregate arrangements were performed. The aggregate 
diameter was varied between 2 and 16 mm. Furthermore, a range of volume fractions between 0.1 and 
0.5 was studied. The nonlinear analyses were based on a 2D lattice approach in which aggregates were 
simplified as monosized cylindrical inclusions. The analysis results were interpreted by means of crack 
length, crack width and change of permeability. The results show that increasing aggregate diameter 
(at equal volume fraction) and decreasing volume fraction (at equal aggregate diameter), increases crack 
width and consequently greatly increases permeability. 

Keywords: Microcracking (B); Interfacial transition zone (B); Transport properties (C); Shrinkage (C); 
Aggregate (D); Lattice modelling. 



1 Introduction 



Drying of cement based composites, such as concrete and mortar, induces cracking if shrinkage of the 
constituents is either internally or externally restrained. For example, non-uniform drying leads to a 
moisture gradient, which results in non- uniform shrinkage of the specimen. Surface regions shrink faster 
than the inner bulk material, which results in surface cracking [EH. E3|. Additionally, shrinkage might be 
restrained by aggregates within the composite Aggregate-restrained shrinkage can lead to micro- 

cracking, which strongly influences the transport properties of the material [l8l [j, [28| . However, the 
evolution of micro-cracks and their dependence on the size and volume fraction of aggregates is not fully 
understood yet. In 28] it was observed that permeability increases with increasing aggregate size at a 
constant aggregate fraction. This result is surprising, since an increase of the aggregate size at a constant 
aggregate fraction is usually accompanied by a decrease in the volume of interfacial transition zones 
(ITZs), which are known to be more porous than the cement paste. One hypothesis is that an increase 
of aggregate diameter at constant aggregate fraction results in an increase of micro-crack width, which 
is closely related to permeability [27| . The objective of this work was to establish whether this size effect 
really occurs. This will undoubtedly enhance the understanding of the link between micro-structure and 
macro-property, in particular the effect of microcracking on mass transport, which is a critical aspect for 
predicting durability and service-life. Shrinkage induced micro-cracking has been investigated recently 
with a lattice model (25[. However, to the authors' knowledge the influence of the size of the aggregates 
on micro-cracking has not been analysed before. 



In the present work, shrinkage induced micro-cracking of concrete was analysed by means of the nonlinear 
finite element method. In concrete, stress-free cracks form by a complex nonlinear fracture process, during 
which energy is dissipated in zones of finite size. There are three main approaches to describe fracture 
process zones within the finite element framework. Continuum approaches describe the evolution of cracks 
as zones of inelastic strains by means of higher-order constitutive models such as integral-type nonlocal 
models [!, 23|. In hybrid approaches, cracks are modelled as displacement discontinuities, which are 



embedded into the continuum description 2(J ]jj 0] ■ Finally, discrete approaches represent the nonlinear 
fracture process by means of the failure of discrete elements, such as trusses and beams 24, ljj. In recent 
years, one type of discrete lattice approach based on the Voronoi tessellation has been shown to be very 
suitable for fracture simulations [5j, |f3j. This lattice approach is robust and computationally efficient. 
With specially designed constitutive laws, the results obtained with this modelling approach are mesh- 
independent. In the present study, this lattice approach was used in combination with a damage-plasticity 
constitutive model, which was designed to result in mesh-independent responses [l(| • 



The present study is based on several simplifications. Shrinkage is represented by an eigenstrain, which 
was uniformly applied to the cement matrix only. This is representative of autogenous shrinkage, but 
does not fully represent transient non-uniform shrinkage due to moisture gradients which may lead to 
more cracking near the surface than in the centre of a concrete element. Furthermore, the only inclusions 
considered were aggregates, which were embedded in a uniform cement paste. Micro-cracking due to other 
inclusions, e.g. unhydrated cement and calcium hydroxide crystals, was not considered. Aggregates were 
assumed to be separated from the cement paste by interfacial transition zones (ITZs) , which were modelled 
to be weaker and more brittle than the cement matrix (l6| . Shrinkage of aggregates, which might have 
a considerable influence on the overall shrinkage of concrete, was not modelled. However, the shrinkage 
strain applied to the cement paste can be interpreted as the difference of uniform shrinkage in cement 
paste and aggregates. The mechanical response of aggregates was assumed to be elastic. Furthermore, 
the study was limited to two-dimensional plane stress analyses, in which aggregates are idealised as 
cylindrical inclusions of constant diameter. 



2 Experimental results 



This section summarises the findings from a recent experimental study by Wong et al. [28[ as background 
to and motivation for the modelling work presented in this paper. Concretes and mortars with a range 
of aggregate contents, w/c ratio, binder type, curing period and preconditioning temperature were tested 
for oxygen diffusivity, oxygen permeability and water sorptivity. Thames Valley gravel (5-12. 7mm) and 
sand (<5mm) were used as coarse and fine aggregates respectively. The main objective of the study was 
to examine the relative influences of ITZ and microcracking induced by oven drying on different transport 
properties, in light of the inconsistencies in previous investigations on ITZ and mass transport. 

Fig. [U reproduces the diffusivity and permeability results for 90-day cured samples that were precondi- 
tioned at 50°C, 10% r.h. Note that this represents only a small portion of the total results presented by 
Wong et al. [28(, but captures the important aspects of their findings. Interested readers are referred 
to the original article for experimental details and in-depth discussion. It was observed that samples 
with low w/c ratio and containing silica fume recorded the lowest transport coefficients for all aggregate 
contents and curing ages, as expected. The transport coefficients also decreased with longer curing age, 
but increased at higher conditioning temperature. For mortars, the transport coefficients decreased sig- 
nificantly with increase in sand content. At very high sand fractions, the ITZs should be overlapping, 
but this did not have any detrimental effect on the transport properties. The trend remains consistent 
regardless of the type of transport property, w/c ratio, binder, curing age and even after severe oven 
drying treatment at 105°C (results shown in [28|). 

The most interesting finding was that concrete has roughly the same diffusivity and sorptivity as its 
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Figure 1: Effect of aggregate volume fraction on a) oxygen diffusivity and b) oxygen permeability, of 
mortars and concretes with different w/c ratios (0.5 & 0.3) and binder type. The samples were sealed 
cured in cling film for the first 90 days and then preconditioned at 50° C, 10% relative humidity for a 
further 90 days until constant mass prior to testing. The dotted circle highlights the anomalous trend 
for concrete samples, which are represented with squares. (M = mortar; C = concrete; SF = containing 
8% silica fume). 



analogous mortar, but significantly higher permeability (1-2 orders of magnitude), despite having about a 
third less ITZ (see Fig. [I]). This striking feature is also consistent across different w/c ratios, binder types, 
curing ages and conditioning regimes, increasing the confidence that it is not caused by experimental errors 
or a series of coincidences. Furthermore, the same sample was used for all three transport tests, and the 
permeability test was performed in between diffusivity and sorptivity, which discounts the possibility that 
the samples were damaged during permeability testing. However, the higher permeability in concretes 
cannot be attributed to the ITZs, because it would have a larger influence on the mortars if this was the 
case. 

Wong et al. [28| suggested that the higher permeability is due to more microcracks forming in concrete 
compared to in the analogous mortar, when both are subjected to the same drying condition. This is 
because aggregates in concrete have larger size distribution, which at the same volume fraction, equates 
to lower surface area to bond with the paste. Thus, any shrinkage (or expansion) in the paste would 
produce larger stresses at the interface and more microcracking compared to the analogous mortar. Since 
diffusivity and sorptivity are less sensitive to microcracking than permeability, these properties are not as 
affected. Microcracks in concretes also have lower specific lengths and tortuosity, because there are fewer 
obstructing aggregates per unit volume to arrest the cracks. Indeed, backscattered electron microscopy 
and image analysis on a few samples found that this may be the case. However, experimental results are 
often affected by many varying parameters that are difficult to isolate and control. Therefore, the present 
numerical study was carried out to confirm the trends observed in the experiments and to increase the 
understanding of the underlying mechanisms. 
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Figure 2: Discretisation: (a) Lattice elements (solid lines) and cross-sections (dashed lines) obtained from 
the Delaunay triangulation and dual Voronoi tessellation, respectively, (b) Degrees of freedom u±, v\, tf>±, 
U2, V2 and <p2 of the lattice element of length h and cross-section width I in the local coordinate system. 
The point C at which the interface model is evaluated is in the center of the polygon facet at a distance 
e from the center of the lattice element. 

3 Modelling approach 

In the present work, shrinkage induced cracking was described by means of a lattice approach Q combined 
with a damage-plasticity constitutive law [l0(. Here, the modelling approach is briefly reviewed. Nodes 
are placed randomly in the specimen constrained by a minimum distance d m , i.e. the smaller <j m is, 
the smaller the average element length (Fig. [5^). Based on these randomly placed nodes, the spatial 
arrangement of lattice elements is determined by a Voronoi tesselation. The cross-sections of the lattice 
elements, which connect the nodes (Voronoi sides), are the edges of the Voronoi polygons (Fig.^k). Each 
node has three degrees of freedom, two translations and one rotation, shown in the local coordinate system 
in Figure [lb- The degrees of freedom u = {m, V\, <j>x, U2, t>2, 02 } T of two nodes of a lattice element are 
related to the displacement discontinuities u c = {u c , v c } T at the mid-point C of the cross-section by 



BU C (1) 



where 

"10 e 1 -e 



B 



-1 -h/2 1 -h/2 



(2) 



The variable h is the element length and e is the eccentricity, defined as the distance between the mid- 
points of the lattice element and the corresponding cross-section, respectively (Figure HJd) . The clement 
stiffness in the local co-ordinate system is 

K c = ^-B T DB (3) 

h 

where I is the length of the cross-section (Figure [2b) and D is the constitutive matrix relating stresses 
to strains. Different material phases are modelled by placing nodes strategically, so that cross-sections of 
lattice elements form the boundary between the phases [29| . 

The constitutive model of the lattice elements is a combination of plasticity and damage mechanics [Io| , 
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normalised strain (e n f t )/(2 G ft ) 



Figure 3: Stress-strain response for fluctuating normal strains for fi — 1 (solid line) and /i = (dashed 
line) . 

which relates stresses to strains. The strains e = (s n ,e s ) T are determined from the displacement jump 
u c = (u D ,u s ) T at mid-point C in Eq. {T]) as 

£ = T (4) 

These strains are related to the nominal stress cr — (<j n7 a s ) T as 

cr = (1 - w) D c (e - £ p - e s ) = (1 - to) a (5) 

where ui is the damage variable, D c is the elastic stiffness, £ p = (£ pn ,£ps) T is the plastic strain, e s = 
(e s , 0) T is the shrinkage eigenstrain and cr = (a n , a s ) T is the effective stress. The elastic stiffness is 

where E and 7 are model parameters controlling both the Young's modulus and Poisson's ratio of the 



material For a regular lattice and plane stress, Poisson's ratio v is 

1-7 
3 + 7 



(7) 



The plasticity part of the model is based on the effective stress er and consists of an elliptic yield surface, an 
associated flow rule, an evolution law for the hardening parameter and loading and unloading conditions. 
A detailed description of the components of the plasticity model is presented in [10( . The initial yield 
surface of the plasticity model is determined by the tensile strength / t , the shear strength s/ t and 
the compressive strength eft- The evolution of the yield surface during hardening is controlled by the 
model parameter /i, which is defined as the ratio of permanent and reversible inelastic displacements. 
The damage part is formulated so that linear stress inelastic displacement laws for pure tension and 
compression are obtained, which are characterised by the fracture energies Gft and Gf c , respectively. 

The constitutive response of the interface model is demonstrated in Fig. [3] by the stress-strain response 
for fluctuating normal strains for /j, — 1 and /1 = 0. The normal strain is increased to point A (A). 
Then the strain is reduced to point B (B ) and again increased to point C (C ). For fi — * 0, a pure 
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(a) (b) 

Figure 4: (a) Geometry of the test specimen of length L. Aggregate size <f> and volume ratio p are varied. 
Specimen length L is determined according to Eq. (|10p . (b) Lattice for numerical analysis for p = 0.3 
and cj> — 16 mm. 

damage-mechanics response is obtained, which is characterised by reversible inelastic strains. For /i = 1, 
on the other hand, a pure plasticity model is obtained. The unloading is elastic and the compressive 
strength is reached sooner than for p = 0. The equivalent crack opening w c is defined by the equivalent 
inelastic displacement, which for the present damage-plasticity model is defined as 

*c = IK|| (8) 

where 

w c = h (e p + u (e - £ p )) (9) 

The inelastic displacement vector w c is composed of a permanent and reversible part, defined as he v and 
hu) (e — £ p ), respectively. 

4 Nonlinear finite element analysis of shrinkage induced micro- 
cracking 

Shrinkage induced micro-cracking was analysed by means of the nonlinear finite element approach de- 
scribed above. The elements representing the cement paste were subjected to an incrementally applied 
uniform shrinkage strain up to e s — 0.5 % (Eq. [5|). This value was chosen for the simulation to represent 
a relatively severe shrinkage of the neat cement paste on first-drying to low humidities. For example, 
Helmuth and Turk [12]] reported first-drying shrinkage values of 0.35-0.70% for neat cement pastes with 
w/c ratios of 0.3-0.6, cured for 6-30 months and dried at 47% r.h. Further drying to 7% r.h. produced 
additional irreversible shrinkage of 0.1%. In another study, Fu et al. [H measured shrinkage eigenstrain 
values of 0.2-1.3% by oven drying neat cement pastes (w/c ratios 0.25-0.70, cured for 6 months) at in- 
creasing temperature up to 110°C. The influence of aggregate volume fraction and aggregate diameter 
was studied. Aggregate volume fractions p = 0.5, 0.3 and 0.1 were modelled. Furthermore, four different 
aggregate diameters {<j> = 16, 8, 4, and 2 mm) were used. The geometry of the specimen analysed for all 
volume fractions and aggregate sizes is shown in Figure Hk. The length L of the specimen in Fig.0^ is 
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Table 1: Model parameters. 



Phase 



E [GPa] 7 ft [MPa] s c G {t [3 /mm 2 } G lc [3 /mm 2 } p 



Cement paste 



40 0.33 6.5 2 20 100 100000 

57.1 0.33 3.25 2 20 50 50000 

100 0.33 - - - - - 



ITZ 

Aggregate 




(10) 



Thus, the smaller the aggregate diameter <f> and the greater the volume ratio p, the smaller is the specimen 
length L. The smallest separation distance between aggregate particles (s a -(/>), i.e. minimum width of the 
cement paste decreases with increasing aggregate volume fraction at constant aggregate size, as would be 
expected in the case of real mortars with increasing aggregate fraction. At constant aggregate fraction, 
the separation between aggregate particles increases with increasing aggregate size, similar to the case 
of a mortar compared to concrete at the same aggregate content. The ratio of the aggregate diameter <p 
and the minimum distance d m for the nodes of the lattice was chosen as 64/3. We will show later that 
varying this ratio has little effect on the simulation results. Thus, the same detail of discretisation was 
applied, independently of the size of aggregates. For a diameter of <fi = 4 mm for instance, the minimum 
distance corresponds to d m = 0.1875 mm. The material parameters for the constitutive model were chosen 
according to Tab. Q] [9]. Aggregates were modelled elastically. Lattice elements crossing the boundary 
between aggregates and cement paste represent the average response of the interfacial transition zones 
and the two adjacent material phases. In all the analyses, the length of the lattice elements is significantly 
greater than the width of the highly non-uniform interfacial transition zones, which is usually in the range 
of 10s of ^m [jjj l. Therefore, the stiffness of these lattice elements was approximated as 



where E a and E m are the Young's modulus of aggregate and mortar respectively. Consequently, the 
Young's modulus of these elements is independent of the element size. The strength of these lattice 
elements is determined by the strength of the ITZ, which is the weakest link. Here, the strength and 
fracture energy ratio of cement paste and ITZ was chosen as 2. This ratio is an approximate value for 
samples with a relatively weak ITZ. For instance, Hsu and Slate [3] observed that the average tensile 
bond strength of the paste-aggregate interface varied from 50-75% of the paste tensile strength, depending 
on the aggregate type, surface roughness and w/c ratio (0.265-0.36). The strength ratio, however, was 
observed to be independent of curing age (3-90 days). 

With p = 0.001, the constitutive model response is close to a perfect elasto-damage response (Fig. [3]), 
which is characterised by the absence of permanent displacements. For an aggregate diameter of 4> — 
16 mm and a volume fraction of p = 0.3, the crack patterns and the deformed mesh are shown in 
Fig. and b. Black lines mark cross-sections of lattice elements, for which the maximum crack width 
is reached at this stage of the analysis. This crack stage is defined as active. On the other hand, grey 
lines represent lattice elements, which reached their maximum crack widths at an earlier stage of the 
analysis. The crack opening in these elements is either reducing, or increasing again after an earlier 
reduction. The simulated widths of the active and passive cracks range between and 70 pm. Some 
bond cracks between the aggregate and paste can be observed, however most cracks appear to originate 
at the aggregate surface and propagate towards the matrix. With the deformed mesh in Fig. it is 
illustrated that the deformations are localised in a few cracks, which connect the aggregates in a regular 
square pattern. Qualitatively very similar crack patterns were obtained for the other aggregate diameters 
and volume fractions. This pattern of cracking seems consistent to that observed by Hsu [l5( in '2D' 
model samples made of sandstone discs arranged in a square grid and filled with paste that is subsequently 




(11) 
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(a) 



(b) 



Figure 5: (a) Crack patterns obtained for p = 0.3 and <fi = 16 mm at a shrinkage strain in the matrix 
of e s — 0.5 %. Black lines indicate active cracks, which increase at this stage of analysis. Grey lines 
represent inactive cracks, which do not increase at this stage of analysis, (b) Deformations (magnified by 
a factor of 50) at the same shrinkage strain. 

subjected to drying shrinkage. Depending on the separation between the aggregates, cracks were seen to 
occur at the interface (i.e. bond cracks), near the shortest distance and diagonally at the largest distance 
between aggregates. 

Length and width of micro-cracks for different volume fractions and aggregate diameters were compared. 
The specific crack length l c is 



where A = L 2 is the area of the specimen, li is the width of the cross-section (Fig [2]) of lattice element i, 
and n c is the number of cracked lattice elements. Furthermore, the average crack width is 



where w C i is the crack width of element i (Eq. [8|) . The influence of the aggregate diameter and volume 
fraction on the specific crack length and average crack width are shown in Figs. and b, respectively. 
The specific crack length in Fig. [(Ji increases with decreasing aggregate size at constant volume fraction. 
An increase of volume fraction at constant aggregate size results in a decrease of the specific crack length. 
Furthermore, the average crack width (Fig. Eh) decreases with decreasing aggregate size at constant 
volume fraction. Again, an increase of the volume fraction at constant aggregate size results in a decrease 
of the crack width. Consequently, a decrease of the aggregate size results in more micro-cracking, which 
is in agreement with experimental observations made in |28j . The crack width is controlled by the spacing 
of the dominant cracks (Fig. which is a function of aggregate diameter and volume fraction. The 
greater the spacing between cracks, the greater is the average crack width. 

Crack width is closely related to transport properties of concrete, in particular in the case of flow under a 
pressure gradient. Assuming that the paste matrix is dense so that flow occurs predominantely through 
the cracks, the permeability k c in the out-of-plane direction due to cracking can be described by the cube 
of the crack width as 




(12) 




(13) 




(14) 
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Figure 6: (a) Specific crack length l c in Eq. (fT2|) versus aggregate diameter <f> for volume fractions of 
p = 0.1,0.3 and 0.5 at a shrinkage strain of e s = 0.5 %. (b) Average crack width w c in Eq. (fT3")) versus 
aggregate diameter <f) for volume fractions of p — 0.1, 0.3 and 0.5 at a shrinkage strain of e s — 0.5 %. 

where £ is a material constant [U El ■ The influence of aggregate diameter and volume fraction on k c is 
shown in Fig.[7K on a log- log scale for £ = 1. Note that what is of interest here is not the actual value of the 
estimated permeability, but the change in permeability caused by varying either the aggregate volume 
fraction or particle size. At a constant aggregate volume fraction, increasing the aggregate diameter 
from 2 to 16 mm caused approximately a 2.5 orders of magnitude increase in permeability. At constant 
aggregate diameter, increasing the aggregate volume fraction from 0.1 to 0.5 produced about 1 order of 
magnitude decrease in permeability. The aggregate diameter influences the permeability strongly, since 
the crack width, which increases with increasing aggregate diameter (Fig. [5b), enters Eq. (TT4")) in its cube. 



In Fig. [5^ many inactive cracks (grey lines) are visible, whereas the displacements are localised in only 
a few zones as illustrated by the deformed specimen in Fig. [5Jd. These inactive cracks have small widths 
(<1 pm), but the evaluation of crack length in Eq. (fl"2|) and crack width in Eq. (I13|) would be strongly 
influenced by this initial cracking. To avoid misinterpretation, only lattice elements with crack-openings 
greater than 1 pm were considered in the evaluation of Eqs. (fT2"f to (|T4"|) . 



5 Discussion 



The edge length L of the specimen (Eq. fTU|) depends on the aggregate diameter <\> and the volume 
fraction p. Thus, the specimen area for which crack length, crack opening and hydraulic conductivity 
were evaluated, differed depending on the micro-structure considered. Comparison of the results for the 
different micro-structures is only valid if the average properties are independent of the specimen area. 
Here it is demonstrated that this is the case. It is assumed that the square crack pattern, shown in 
Fig- Eh in the form of cracks connecting the aggregate particles with a spacing s a , is independent of the 
number of aggregate particles considered. Thus, the number of dominant cracks in one direction of the 
specimen is n c — L/s a . For instance, for the crack patterns shown in Fig. [5Jd, n c — 2, whereas for a 
specimen of double the size it would be 4, since s a is related to the material structure and is independent 
of the specimen length L. With this assumption, the specific crack length l c in Eq. (fT2|) . which is defined 
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Figure 7: (a) Crack permeability k c (Eq. [T4|) versus aggregate diameter </> for volume fractions of p = 
0.1,0.3 and 0.5 at a shrinkage strain of e s = 0.5 %. (b) Minimum spacing of lattice nodes versus 
permeability fc c for p = 0.3 and = 4 mm. 

as the total length of cracks divided by the specimen area, is 



_ 2 2 ( 1 2<fi 



L — - 1 — 

2n c (L - n c (j)) _ s a \ s a J _ 2 f l 20 



■5a \ 



(15) 



where L — n c <\> is the length of one crack. The specific crack length l c depends only on the parameters 
s a and <p, but not on the specimen length L. An increase of L, at constant <f> and s a , does not change l c 
if the influence of the specimen boundaries, which is studied later, is negligible. Equivalent observations 
can be made for the average crack width w c and k c . 

The damage-plasticity constitutive model used for the finite element analyses was designed to result 
in mechanical responses, which are independent of the size of the lattice elements flpj ] . In the present 
study, a mesh-size independent description of the crack opening is important, since the analysis for 
varying aggregate diameters for the aggregate arrangement in Fig. was performed with the same 
diameter mesh size ratio of 64/3. Consequently, the smaller the aggregate diameter, the smaller is the 
length of the lattice elements used. To investigate a possible influence of element size on crack-openings, 
two additional analyses with 4> to d m ratios of 32/3 and 128/3 were performed for a volume fraction 
p = 0.3 and a constant aggregate diameter <fi = 4 mm. The influence of element size on permeability 
(Eq. fT4"|) with £ = 1 is shown in Fig. 03 on log-log scale. The evaluation of permeability is almost 
independent of the element size. The change in the simulated permeability due to element size is small 
and negligible compared to the change in permeability due to varying aggregate size and volume fraction 
(Fig-IZtO- As mentioned above, the results of the analyses were evaluated for lattice elements with crack 
openings greater than 1 /im. This threshold was applied to avoid misinterpretation due to distributed 
initial cracking. The influence of the crack opening threshold on crack length (Eq. I12p and permeability 
(Eq.[T4|) with £ = 1 is shown in Fig.EK and b, respectively. The specific crack length is strongly influenced 
by the crack threshold. On the other hand, crack permeability is almost independent of the threshold, 
since lattice elements with small crack openings contribute little to the overall crack permeability. 

The meso-structure used in this study is a strong idealisation of concrete and results in regular square crack 
patterns (Fig. Eb), which differ considerably from the random crack patterns observed in experiments 
[28j. Therefore, additional analyses were carried out to investigate the influence of random aggregate 
arrangements for a volume fraction of p = 0.3 and two aggregate diameters of (f> — 4 and 8 mm. The 
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Figure 8: Influence of crack opening threshold for p = 0.3 and <j> — 4 mm. (a) Specific crack length versus 
crack opening threshold, (b) Crack permeability k c (Eq. 1 14|) versus crack opening threshold. 

specimen length was chosen as L = 50 mm. Three analyses for each aggregate diameter were carried out 
to determine the mean of the random results. The same material properties as for the regular analyses 
were used. The crack patterns for two random analyses for aggregate diameters = 4 and 8 mm are 
shown in Fig. [9^ and b, respectively. The specific crack length l c , the average crack width w c and the 
change of permeability k c obtained by means of the random analyses are presented in Fig. [6^ and b, 
and Fig [TJi (triangular symbols) together with the results obtained from the analysis with regular ag- 
gregate arrangements. Although the resulting crack pattern differs considerably, as would be expected, 
the specific crack length, average crack width and permeability results for random and regular aggre- 
gate arrangements are in good agreement. Consequently, the regular arrangement seems appropriate to 
investigate the influence of aggregate diameter on shrinkage induced cracking. 

It may be useful to relate the findings of this study to that of the experimental results by Wong et al. 
[28l ]. A direct comparison is not possible because of the approximations used in this study, although 
certain trends seem to be captured by our simulations. Using image analysis, Wong et al. [28( found that 
the concretes have more microcracks, and that the microcracks have lower specific lengths, compared 
to analogous mortars subjected to the same drying conditions. Fig. [TO] shows crack patterns observed 
on a mortar containing 40% vol. sand (0.5 w/c ratio, 90-day cured) after drying at 50°C, 10% r.h. 
using backscattered electron microscopy. Examination of samples dried at 105°C found more severe 
cracking, and significantly higher permeability [28}. The microcracks have widths of about 0.5-10/im, 
although finer cracks were seen at higher magnifications. The crack pattern bears some resemblance to 
the simulations shown in Fig. [9] such as its random orientation, the partial bond cracks appearing around 
some aggregate particles and the matrix cracks that propagate through the paste and very often bridging 
several aggregate particles. However, the simulations did not produce cracks that propagate through the 
aggregates. Also, the simulated cracking appears to be more severe than that observed on the mortars 
dried at 50°, 10% r.h. A possible reason for this is a mismatch in the model parameters used for the 
simulation. 

Wong et al. [28| found that increasing aggregate fraction of the mortar from 10% to 50% resulted in a 
decrease in measured permeability by about an order of magnitude (Fig. [TJs). Given the same drying 
treatment, concretes (with max. aggregate size 12.7 mm) have significantly higher permeability, about 
1.5 orders of magnitude, than mortars (max. aggregate size 5 mm) at the same aggregate content. 
Interestingly, this change in permeability due to the aggregate content and size effect is close to the 
simulated values in Fig. [7J This obviously may be a fortuitous result since several important features 
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Figure 9: Crack patterns for p = 0.3 for (a) (f> = 4 mm and (b) (j) — 8 mm at a shrinkage strain in the 
matrix of e s = 0.5 %. Black lines indicate active cracks, which increase at this stage of analysis. Grey 
lines represent inactive cracks, which do not increase at this stage of analysis. 




Figure 10: Segmented BSE images of a 90-day mortar (w/c 0.5, 40% vol. sand) that was conditioned 
at 50°C, 10%r.h. showing the microcracks, pores, aggregates and cement paste. The image on the left 
is a composte image produced by aligning 3x3 images captured at 100 x, to give a larger field of view 
(3100x400/im). The image on the right was obtained at 200x. Pores with aspect ration < 2 were filtered 
to highlight the microcracks. 
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were not captured by the model, for instance, the contribution of capillary pores, the effect of tortuosity 
and connectivity, shape and size distribution. Some may have an opposite effect or a minor contribution 
to transport properties, but clearly more work needs to be carried out. Nevertheless, the results of this 
study show that the aggregate content and particle size have a significant influence on the formation of 
microcracking induced by aggregate restrained shrinkage, which in turn affects the permeability of the 
composite. 



6 Conclusions 

In the present work the influence of aggregate size and volume fraction on shrinkage induced micro- 
cracking was studied numerically by means of the nonlinear finite element method. The work resulted in 
the following conclusions: 

• The length of micro-cracks decreases with increasing aggregate diameter and increasing volume 
fraction. 

• The average crack width increases with increasing aggregate diameter and decreasing volume frac- 
tion. 

• The permeability, which is related to the cube of the crack width, increases with increasing aggregate 
diameter and decreasing volume fraction 



Thus, the aggregate diameter and volume fraction influence the formation of microcracks induced by 
aggregate restrained shrinkage, which is in agreement with experimental observations in [281 ] (Sec. [5]). 
However, the representation of aggregates as circular inclusions and the use of one aggregate size over- 
simplifies the micro-structure of concrete, which does not allow a direct comparison with experimental 
results. In future work, the present modelling approach will be extended to 3D and the use of realistic 
aggregate grading curves will be included by applying multi-scale analysis approaches. 



7 Acknowledgements 

The simulations were performed with the object-oriented finite element package OOFEM 21 , 22| extended 
by the present authors. HSW and NRB acknowledge the financial support from the Engineering and 
Physical Sciences Research Council, UK. 



References 

[1] CM. Aldea, M. Ghandehari, S.P. Shah, and A. Karr. Estimation of water flow through cracked 
concrete under load. ACI Materials Journal, 97(5):567-575, 2000. 

[2] Z. P. Bazant and W. J. Raftshol. Effect of cracking in drying and shrinkage specimens. Cement and 
Concrete Research, 12(2):209-226, 1982. 

[3] Z.P. Bazant and M. Jirasek. Nonlocal Integral Formulations of Plasticity and Damage: Survey of 
Progress. Journal of Engineering Mechanics, 128:1119, 2002. 

[4] J. Bisschop and J.G.M. van Mier. Effect of aggregates on drying shrinkage microcracking in cement- 
based composites. Materials and Structures, 35(8):453-461, 2002. 



13 



[5] J. E. Bolander and S. Saito. Fracture analysis using spring networks with random geometry. Engi- 
neering Fracture Mechanics, 61:569-591, 1998. 

[6] J. E. Bolander and N. Sukumar. Irregular lattice model for quasistatic crack propagation. Physical 
Review B, 71, 2005. 

[7] G. T. Camacho and M. Ortiz. Computational modelling of impact damage in brittle materials. 
International Journal of Solids and Structures, 33:2899-2938, 1996. 

[8] Y. Fu, P. Gu, P. Xie, and J. J. Beaudoin. Development of eigenstress due to drying shrinkage in 
hardened portland cement pastes: Thermomechanical analysis. Cement and Concrete Research, 24 
(6):1085-1091, 1994. 

[9] P. Grassl. Meso-scale analysis of concrete splitting tests. Internal report GUCE2008PG02, Depart- 
ment of Civil Enginering, University of Glasgow, Glasgow, United Kingdom, 2008. 

[10] P. Grassl and R. Rempling. A damage-plasticity interface approach to the meso-scale modelling of 
concrete subjected to cyclic compressive loading. Engineering Fracture Mechanics, 75:4804-4818, 
2008. 

[11] D. V. Griffiths and G. G. W. Mustoc. Modelling of elastic continua using a grillage of structural 
elements based on discrete element concepts. International Journal for Numerical Methods in Engi- 
neering, 50:1759-1775, 2001. 

[12] R. J. Hclmuth and D. H. Turk. The reversible and irreversible drying shrinkage of hardened portland 
cement tricalcium silicate pastes. Journal of Research and Development Laboratory of Portland 
Cement Association, 9:8-21, 1967. 

[13] H.J. Herrmann, A. Hansen, and S. Roux. Fracture of disordered, elastic lattices in two dimensions. 
Physical Review B, 39(l):637-648, 1989. 

[14] D. W. Hobbs. Influence of aggregate restraint on the shrinkage of concrete. ACI Journal Proceedings, 
71(9), 1974. 

[15] T. T. C. Hsu. Mathematical analysis of shrinkage stresses in a model of hardened concrete. A CI 
Journal, 60(3):371-390, 1963. 

[16] T. T. C. Hsu and F. O. Slate. Tensile bond strength between aggregate and cement paste or mortar. 
A CI Journal Proceedings, 60(4), 1963. 

[17] CL Hwang and JF Young. Drying Shrinkage of Portland Cement Pastes I. Microcracking During 
Drying. Cem. Conor. Res, 14:585-594, 1984. 

[18] A.D. Jensen and S. Chatterji. State of the Art Report on Microcracking and Lifetime of Concrete- 
Part I. Materials and Structures, 29(l):3-8, 1996. 

[19] M. Jirasek. Comparative study on finite elements with embedded cracks. Computer Methods in 
Applied Mechanics and Engineering, 188:307-330, 2000. 

[20] N. Moes and T. Belytschko. Extended finite element method for cohesive crack growth. Engineering 
Fracture Mechanics, 69:813-833, 2002. 

[21] B. Patzak. Object oriented finite element modeling. Acta Polytechnica, 39:99-113, 1999. 

[22] B. Patzak and Z. Bittnar. Design of object oriented finite element code. Advances in Engineering 
Software, 32:759-767, 2001. 

[23] G. Pijaudier-Cabot and Z. P. Bazant. Nonlocal damage theory. Journal of Engineering Mechanics, 
ASCE, 113:1512-1533, 1987. 



14 



[24] E. Schlangen and J. G. M. van Mier. Simple lattice model for numerical simulation of fracture of 
concrete materials and structures. Materials and Structures, 25:534-542, 1992. 

[25] E. Schlangen, EAB Koenders, and K. van Breugcl. Influence of internal dilation on the fracture 
behaviour of multi-phase materials. Engineering Fracture Mechanics, 74(l-2):18-33, 2007. 

[26] K.L. Scrivener, A.K. Crumbie, and P. Laugesen. The interfacial transition zone (ITZ) between 
cement paste and aggregate in concrete. Interface Science, 12(4):411-421, 2004. 

[27] P. A. Witherspoon, J. S. Y. Wang, K. Iawai, and J. E. Galw. Validity of cubic law for fluid flow in 
a deformable rock fracture. Water Resour. Res, 16(6):1016-T024, 1980. 

[28] H. S. Wong, M. Zobel, N. R. Buenfeld, and R. W. Zimmerman. The influence of the interfacial 
transition zone and microcracking on the diffusiviy, permeability and sorptivity of cement-based 
materials after drying. Magazine of Concrete Research, 2008. In press. 

[29] M. Yip, J. Mohle, and J. E. Bolander. Automated Modeling of Three-Dimcnsional Structural Compo- 
nents Using Irregular Lattices. Computer- Aided Civil and Infrastructure Engineering, 20(6):393-407, 
2005. 



15 



